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ABSTRACT 

The next generation of gamma-ray telescopes may be able to observe gamma-ray 
blazars at high redshift, possibly out to the epoch of reionization. The spectrum of 
such sources should exhibit an absorption edge due to pair-production against UV 
photons along the line of sight. One expects a sharp drop in the number density of 
UV photons at the Lyman edge ei- This implies that the universe becomes transparent 
after gamma-ray photons redshift below E ~ (m e c 2 ) 2 /eL ~ 18GeV. Thus, there is only 
a limited redshift interval over which GeV photons can pair produce. This implies that 
any observed absorption will probe radiation fields in the very early universe, regardless 
of the subsequent star formation history of the universe. Furthermore, measurements 
of differential absorption between blazars at different redshifts can cleanly isolate the 
opacity due to UV emissivity at high redshift. An observable absorption edge should 
be present for most reasonable radiation fields with sufficient energy to reionize the 
universe. Lya photons may provide an important component of the pair-production 
opacity. Observations of a number of blazars at different redshifts will thus allow us to 
probe the rise in comoving UV emissivity with time. 

1. Introduction 

Our knowledge of UV radiation fields and energy injection into the IGM at z > 5 is fairly 
tenuous. There are two main constraints: observations of the integrated background light (Madau 
&; Pozzetti 2000, Bernstein et al 1999), and the fact that no Gunn- Peterson trough is observed in 
the spectra of the highest-redshift quasar to date (Fan et al 2000), implying that the universe must 
be reionized by z = 5.8. The upcoming Next Generation Space Telescope (NGST) will be able to 
image high-redshift star clusters or AGNs in rest frame UV continuum emission (Haiman & Loeb 
1997,1998), and their redshifts may be obtained via Ha observations (Oh 1999). Nonetheless, the 
redshift-binned number counts will be fairly sparse, and one is unlikely to probe sufficiently far 
down the luminosity function to get a good measure of the comoving emissivity as a function of 
redshift. 
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Observations of gamma-ray blazars ("grazars") probe extragalactic IR and UV radiation fields, 
by observing the pair production opacity to 7 rays at the high energy end (Gould & Schreder 1967, 
Stecker, De Jager & Salamon 1992, Madau & Phinney 1996, Primack et al 1999). All theoretical 
models have confined their predictions to low redshift grazars, with the exception of Salamon & 
Stecker (1998), who computed the 7-ray opacity up to z=3. They concluded that because the stellar 
emissivity peaks between z=l and z=2, the 7-ray opacity shows little increase at high redshift, and 
thus is not dependent on the initial epoch of galaxy formation. 

To date, EGRET has detected 66 gamma-ray loud blazars (Hartman et al 1999), out to redshifts 
z > 2. The next generation of gamma-ray telescopes (GLAST, CELESTE, STACEE, MAGIC, 
HESS, VERITAS, and Milagro) should greatly enlarge this sample. If the low redshift correlation 
between black hole mass and bulge mass (Magorrian et al 1998) continues to high redshift, then 
it is possible that high-redshift halos could host mini-quasars (Haiman & Loeb 1998, Haehnelt, 
Natarajan & Rees 1998), which should be detectable in rest frame UV emission by NGST and X- 
ray emission by Chandra (Haiman & Loeb 1998, 1999) in the redshift range z ~ 5 — 15. This raises 
the exciting possibility that grazars could be detected at similarly high redshifts. It is worth noting 
that EGRET has detected ~ 56 sources at high Galactic latitudes b > 10° (Mukherjee, Grenier 
& Thompson 1997), with no known counterparts at other wavelengths. Their spatial distribution 
and log N- log S plot can be well fit by a Galactic component plus an isotropic, extragalactic 
contribution. Some of these may well be unidentified high-redshift blazars. 

In this paper, I point out that if grazars are detected at high redshifts z > 3, the pair production 
opacity to gamma ray photons can be used to probe the comoving emissivity longward of the Lyman 
break at these extremely high redshifts, independent of the star formation rate at lower redshifts. 
Due to the small escape fraction of ionizing photons f esc < 5% from host galaxies, as well as the 
high photoelectric opacity of the IGM at these wavelengths, the comoving number density of UV 
photons exhibits a sharp drop at the Lyman edge at all redshifts. Thus, there is only a limited 
pathlength over which a gamma-ray photon can pair produce against UV photons, before it redshifts 
to energies which require UV photons above the Lyman edge for pair production to take place. For 
z < Zbreak, the universe becomes optically thin to the gamma-ray photon. Thus, the detection of 
an absorption edge in a high-redshift grazar places an immediate constraint on the mean radiation 
field over the redshifts zbreak < z < z s . Furthermore, measurement of the different absorption at 
a given observed energy between blazars at redshifts z±,Z2 places an immediate constraint on the 
radiation field in the redshift range z\ < z < zi- Detection of grazars at a number of redshifts 
would then enable one to probe the UV emissivity history of the universe. 

In all numerical estimates, we assume a background cosmology given by the 'concordance' 
values of Ostriker & Steinhardt (1995): (Q m , Q A , Q b , h, a 8h -i , n) = (0.35, 0.65, 0.04, 0.65, 0.87, 0.96). 
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2. Gamma-Ray Blazars 

A detailed study of the detectability of high-redshift blazars is beyond the scope of this paper. 
In this section, I merely show that it is plausible that GLAST will be able to detect high redshift 
blazars. 

With a point source sensitivity of S(E > lOOMeV) ~2x 10 _7 photons s _1 cm~ 2 , EGRET has 
detected ~ 66 high-redshift blazars out to z > 2 (Hartman et al 1999). The associated gamma- 
ray luminosities correspond to L 7 = 10 46 — 10 49 ergs _1 , and typically dominate the bolometric 
luminosity of the AGN, with L^/Lb ~ 1 — 1000. The upcoming gamma-ray telescope GLAST (see 
http:/ /glast. gsfc.nasa.gov) will be 2 orders of magnitude more sensitive, with a detection threshold 
of S(E > lOOMeV) k 2 x 10~ 9 photons s _1 cm -2 for a 5 a detection with a 50 hour integration, 
and S(E > lGeV) ~ 10~ 10 photons s _1 cm -2 (these thresholds correspond to the same detection 
limit for a L„ oc v~ a source spectrum where a = 1). Goals for GLAST include a broad energy 
coverage from lOMeV— > 300 GeV, with a spectral resolution of ~ 2% in the > 10 GeV range, a 
field of view of > 3 sr, and a source location determination accuracy of 30 arcsec- 5 arcmin. During 
its lifetime, it will perform an all-sky survey similar to that conducted by EGRET. Sources of the 
same or somewhat fainter luminosity as those detected by EGRET may be seen by GLAST out to 
extremely high redshifts, z ~ 10. 

Will such luminous sources will be present at high redshift? If the AGN is assumed to emit all 
its energy at gamma-ray wavelengths at the Eddington luminosity, the inferred black hole mass is 
extremely high, Mt,h = lO lo M (L e( i ( i/lO 48 ergs _1 ). However, there are two reasons why gamma ray 
sources of high apparent luminosity do not require such massive black holes: (i) even if all photons 
are radiated isotropically, if most of the radiation emerges at high energies (as appears to be the case 
in gamma-ray blazars), then Klein-Nishina effects must be taken into account (Dermer & Gehrels 

1995). The inferred black hole mass, given by M^ N > 2 l ^gx^o^ergs- 1 ^l+l)^ 11 ^ 61 ^ + z )] ( wnere 
F(ei,e u ) is the observed flux between the lower and upper bandpass limits ei,e u ) typically drops 
by 3 orders of magnitude, so a lO^ergs" 1 source only requires a 1O 7 M black hole, (ii) There is 
strong evidence for relativistic beaming in blazars (e.g., through the observation of superluminal 
jets (von Montigny et al 1995)). In fact, if blazars were not beamed, we would not be able to see 
any gamma-rays from them due to the high pair production opacity at the source; beaming reduces 
the luminosity /radius ratio by a factor S p+1 , allowing photons to escape (Maraschi et al 1992). 
Here, if a is the spectral index of the source, then p = 3 + a for a moving sphere in the Synchrotron 
Self-Compton (SSC) model, while p = 4 + 2a in the External Radiation Compton (ERC) model; 
5 = [7(1 — PcosO}^ 1 is the relativistic Doppler factor, and 7 is the Lorentz factor. If C is the initial 
intrinsic luminosity of the jet in gamma-ray emission, beaming boosts the observed luminosity of 
the jet to L = 5 P C. For 6 ~ 0°, then 5 ~ 27, and the observed luminosity is amplified by a factor of 
thousands. The strong relativistic beaming reduces the fraction of sources which are visible, since 
they can only be seen when viewed along the jet axis. For instance, for 7 = 6, and a = 1, in the SSC 
model the observed luminosity is reduced by an order of magnitude from its maximum if the jet is 
pointing 8.5° from our line of sight, and two orders of magnitude if the jet is pointing 14.2° from 
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our line of sight. Note that the black hole masses derived for a number of low redshift blazars from 
variability timescale and transparency arguments lie in the range 10 7 - 1O 8 M (Cheng et al 1999, 
Hartman et al 1996, Becker &; Kafatos 1995, Romero et al 2000). The luminosity of these blazars is 
so high they could be seen at high redshift with GLAST, and black hole masses of 10 — 10 8 M Q are 
reasonably abundant at high redshift in certain models of AGN formation (Haiman & Loeb 1998). 
Note that for a set of blazars of constant intrinsic luminosity C and comoving number density, the 
redshift distribution of detected sources in a flux-limited survey flattens considerably and extends 
to higher redshifts as the Lorentz factor increases (Dermer & Gehrels 1995). 

At present, the modelling of even the low-redshift population of gamma-ray blazars is a matter 
of considerable debate. Models which attempt to account for the unresolved gamma ray background 
with faint blazars either extrapolate the observed 7-ray luminosity function obtained with EGRET 
(Chiang & Mukherjee 1998) or use an assumed conversion between the observed radio loud AGN 
luminosity function and the blazar luminosity function (Stecker Sz Salamon 1996). In this paper, 
I use a highly simplified model to estimate the detectability of high-redshift blazars. I assume 
that the intrinsic luminosity C of the jet in gamma-ray emission (prior to beaming) scales with 
the luminosity of the accretion disk C = fL^ isk , where the optical B-band luminosity is taken to 
be an accurate reflection of Ldisk (i n particular, assuming the median quasar spectrum of Elvis 
et al (1994), a 1 Mq black hole shining at the Eddington luminosity has a B-band luminosity of 
5.7 x 10 3 Lb,©). Such a jet-disk correlation is observed in the ratio of observed radio (i.e., after 
beaming) to optical luminosities (Falcke , Malkan & Biermann 1995). I also assume relativistic 
beaming with L = 5 P C, where p = 3 + a (and a = 1, the average spectral index observed in the 
EGRET blazars). The change in the observed luminosity function due to the effects of relativistic 
beaming is given by (Urry & Padovani 1995): 

$ obs = J dCP(L\C)<S> mtr (£) (1) 

where the probability of observing luminosity L given the intrinsic luminosity C is given by: 

P{L\C) = P{S)^r = J—cVPL-tP+V/P (2) 
dL /ryo 

I use the fit to the observed B-band luminosity function (J)(Lb,z) from Pei (1995). I adjust the 
relation C = fL^ isk and the Lorentz factor 7 (note that since the distribution of blazar Lorentz 
factors is unknown, I assume they all have the same Lorentz factor) to fit the number of blazars 
detected by EGRET and their redshift distribution. I ignore the effects of blazar flaring, which 
increases the luminosity by some factor A (where typically A ~ 5) some fraction £ of the time 
(where typically £ ~ 0.03), which results in a second term ^<f>(L/A), since these two degrees of 
freedom, the normalization (chosen by selecting 7 and thus the beaming angle) and luminosity 
boost (choosen by the combination L = S p fL^ isk ), are already present in our model. I find that 
£13 = 3.2 x 10~ 2 Lg\ 3 (where L i3 = (L/1O 13 L ) and 7 = 6 provides a good fit (see Fig 1, top 
panel). Note that since £7,13 < (2j) p L°^ 9 13 ~ 600L^ 9 13 (corresponding to 9 = 0), these relations 
result in a SED which is in reasonable agreement with the observed SEDs of gamma-ray blazars 
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(see Fig 1 in Ghisellini et al 1998), where L 1 /Lb ~ 10 — 100 typically, although it can range from 
1 — 1000. Furthermore, the Lorentz factor 7 ~ 6 is in reasonable agreement with the relativistic 
Doppler factor 8 < 2j from models which take into account the SED, time variability, and gamma- 
ray transparency of blazars; the derived 8 ~ 10 — 20 (Ghisellini et al 1998). The somewhat lower 
Doppler factors I have adopted conservatively underestimate the number of high-redshift blazars. 

I then extrapolate this model to high redshift by applying it to the Press-Schechter based 
model of Haiman & Loeb (1998) for quasars. In this model, which is calibrated to the observed 
luminosity function of Pei (1995) at lower redshifts, each halo with T v ; ir > 10 4 K hosts a black hole 
with mass = 10 _3 - 2 M/j a ; o which shines at the Eddington luminosity for t Q ~ 10 6 years. The 
result is shown in the bottom panel of Fig 1, which shows that gamma-ray blazars may be detected 
out to z ~ 7. By the time GLAST is launched, a large database of quasars with known redshifts 
will be available (e.g. from the SLOAN digital sky survey, SDSS (York et al 2000)), and many high- 
redshift blazars can be selected simply by identifying their optical counterparts. I emphasize once 
again that this highly simplified model is only intended to serve as a plausibility argument. The 
main point is that while the Press-Schechter formalism predicts that massive halos Mhaio > 1O 11 M 
expected to host supermassive black holes of the requisite luminosity become exponentially rare at 
high redshift, processes which increase the luminosity of a lower luminosity population (beaming 
L = 5 P C, flaring L = AL) create a power-law tail of bright sources. I have neglected a tail to 
the distribution of Lorentz factors, or flaring, which could further flatten the redshift distribution 
of detectable sources, increasing the maximum redshift out to which sources can be seen. Finally, 
gravitational lensing could bring otherwise undetectable sources into view, although the low optical 
depths for strong lensing (e.g. r ~ 6 x 10~ 3 for z s = 7, Porciani & Madau 2000) imply that this 
should only have a small impact on number counts. 



The pair-production optical depth for a photon observed at energy E Q and emitted from a 
source at redshift z s is given by (e.g., Madau & Phinney 1996): 

s f Zs , di r 1 „ „ w „ x r° , , , ,„ 



where E = (1 + z)E Q , and is given by the criterion that pair production can take place, which 
requires that Ee(l — cos#) > 2(m e c 2 ) 2 . For a given energy E, the pair production cross-section 
a(E,e,6) rises sharply from the threshold energy e t h, reaches a peak of 0.26ar at e = 2e t h, and 
finally falls off as e _1 for e 3> e^. To calculate the optical depth of the universe to gamma-ray 
photons emitted at high redshift, we need to know the number density of photons as a function of 
energy and redshift. Before doing so in detail, I make some simple estimates. 

One can perform a simple order of magnitude estimate to show that the minimal comoving 
emissivity required to reionize the universe implies a high pair-production opacity for E > E^, but 



3. Calculating pair production opacity 
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a low opacity for E < E t h, where we define the threshold energy in the rest frame of the source 
for pair production against UV photons at the Lyman edge as E t h ~ (m e c 2 ) 2 /e£ ~ 18GeV. This 
is because the number density of UV photons plummets at wavelengths shortward of the Lyman 
limit; the pair production opacity is thus dominated by soft photons e < = 13.6eV. For the 
universe to be reionized, a minimum of 1 ionizing photon per baryon must be emitted. Thus, the 
comoving number density of ionizing photons is rij(E > El) > n&(l + n rec ) ~ 1CT 7 (1 + n rec ) cm~ 3 , 
where (1 + n rec ) is the average number of times each baryon recombines during the reionization 
epoch. The pair production cross section peaks when Ee ~ (m e c 2 ) 2 , with a value of a pp ~ 0.26o"t- 
Thus, across a Hubble volume, the optical depth for photons with E < E t h is r ~ n^ tProper (E > 
El)c pp Ih ~ 5x 10 _3 (l+z/10) 1 ' 5 (l+n rec ), which is undetectably small. What is the pair production 
opacity for photons with E > E t h? We can scale the number density of photons longward of 
the Lyman limit with respect to the number of ionizing photons produced. There are 3 factors 
involved: (i) f esc . Observations of our own Galaxy (Dove, Shull, & Ferrara 2000, Bland-Hawthorn 
& Maloney 1999) and local starbursts (Leitherer et al 1995) find the escape fraction of ionizing 
photons from galaxies into the IGM is of order f esc ~ 3 — 6%, and radiative transfer calculations 
find that the escape fraction should decrease strongly with redshift (Woods & Loeb 1999, Ricotti 
& Shull 2000) . The actual comoving number density of ionizing photons produced in starbursts is 
/esc n b(l + n rec), where f esc ~ 5% is the escape fraction of ionizing photons from the halo, (ii) /break- 
Photoelectric absorption in stellar atmospheres produces a sharp drop in the flux at the Lyman 
edge. For a Salpeter IMF with metallicity ~ 10~ 2 Z@, there are f break ~ 5 times as many photons 
emitted longward of the Lyman limit as there are shortward of it, integrated over the history of 
the starburst. (iii) /opacity, due to intergalactic absorption. The energy density U 1 ~ ^jrJv where 
J v ~ tuKnpf where X m f p is the mean free path of an photon of frequency v. While the universe 
is optically thin to photons longward of the Lyman limit, it is optically thick to ionizing photons, 
which have a much shorter mean free path. The energy density of ionizing photons is lower by 
a factor /opacity = (A m / P (E>E^)) ^ 10 due to their high absorption rate. The mean free path of 
ionizing photons decreases strongly with redshift and for z > 2 the radiation field is largely local, 
due to the large increase in the number of adsorbers (e.g., Madau, Haardt & Rees 1999). Thus, 

We haVe £S>Sj ~ feslfbreakf opacity ~ 10 3 (fe ) ~* (^) . This implies that the pair 

production opacity longward of 6l over a Hubble length is much greater, r ~ 5(l+z/10) L5 (l+n rec ). 
A similar estimate may be obtained by normalizing to the observed metallicity of the IGM at z ~ 3, 
Z ~ 1O~ 2 Z ; this can be shown to correspond to ~ 10 ionizing photons per IGM baryon (Miralda- 
Escude & Rees 1997), or a comoving number density of ~ 50n& for photons longward of the Lyman 

break, compared to the previous estimate of ~ 100 ^o^t) ( ^" 5 eafc ) n b- 

The expectation of a sharp drop in the intensity of the ambient intergalactic radiation field at 
the Lyman edge is fairly robust. A recent composite spectrum of 29 Lyman break galaxies (LBGs) 
exhibit an observed flux ratio L(1500)/L(900)=4.6 ± 1.0, which is consistent with no internal pho- 
toelectric absorption, i.e. f esc ~ 1 (Steidel, Pettini h Adelberger 2000). Note that the authors 
themselves stress this result should be treated as preliminary; the result could be due to a large 
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number of uncertainties or selection effects, among them the fact that these galaxies were selected 
from the bluest quartile of LBGs. Even so, if this is typical of all high redshift galaxies, the observed 
Lyman break (likely due to absorption by stellar atmospheres) and expected processing by inter- 
galactic absorption at high redshift imply a Lyman edge ^ 7 [^^| > 50. Similar considerations of 
internal and intergalactic photoelectric absorption apply if quasars dominate the ionizing radiation 
field. 

As a photon redshifts, it has to interact with higher energy photons to pair produce. A photon 
is able to pair produce until it redshifts below E t h, i.e., the universe is optically thin to a photon 
once it redshifts below (1+zth) ~ ^-(l+Zs); where E s is the original energy of the photon at source 
redshift z s . Thus, there is a fairly well-defined pathlength 51 = ((1 + z s )~ 3 / 2 — (l + z^) -3 / 2 ) over 
which a photon may pair produce. This has the fortunate consequence that low redshift photons 
do not affect photons emitted near E t h', thus, near the threshold energies one is always probing 
the average radiation field at redshifts comparable to the source redshift. The flux decrement at 
a given frequency measures the average number density of photons redward of the Lyman break 
over the associated redshift interval, n 7 ~ T/(Ala pp ). Gamma ray absorption measurements thus 
provide a fairly clean measurement of radiation fields at high redshift which are uncomplicated 
by radiative transfer effects since (apart from unimportant H2 opacity effects) the universe is 
optically thin to photons redward of the Lyman limit. In particular, photons longward of the 
Lyman limit establish a homogeneous, isotropic radiation field early in the history of the universe. 
As the pathlengths for pair production opacity to become significant are typically of order a Hubble 
volume, opacity fluctuations due to Poisson fluctuations or source clustering are insignificant. The 
measured background radiation field may be compared directly against the expected background 
from measurements of the source luminosity function by direct imaging of high redshift sources in 
mid-IR with NGST, where fast photometric redshifts may be obtained using the Gunn-Peterson 
break. A comparison of the two should in principle allow one to check the completeness of a survey 
at NGST flux limits. 

Let us now use a specific model of high-redshift star formation to calculation the abundance 
of UV photons at high redshift. The solution of the cosmological radiative transfer equation yields 
the mean specific intensity of the radiation background at the observed frequency v , as seen by an 
observer at redshift z Q as (Peebles 1993): 

JK ,*) = ^**<^(,,, )e -^~ » (4) 

where v = v (\ + z)/(l + z ), and r e ft(^ , z , z) is the effective optical depth of the IGM to radiation 
emitted at z and observed at z Q at frequency v Q . Redward of the Lyman edge, radiative transfer is 
particularly simple as the universe is optically thin, and only the redshifting of photons is important 
(there is one caveat to this statement: the optical depth of the IGM in the Lyman resonance lines 
prior to reionization is very large; thus whenever a photon redshifts into a Ly/3 or higher order 
Lyman resonance, it is reprocessed into a Lya and Balmer or lower order line photon (see Haiman, 
Rees & Loeb 1997)). However, this merely causes a modulation in the spectrum in the 11.2-13.6 
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cV range, redistributing photons to the Lya and Balmer wavelengths. The large typical energy 
intervals of target photons (see Fig 3, top panel) implies that this redistribution causes the pair- 
production opacity to remain the same or increase. I therefore ignore this complication). The 
number density of photons in an energy interval is then given 

b y 7ET = T^ J »- 1 model the star 
formation history fi* of the high redshift universe with the semi-analytic models of Haiman & Loeb 

(1997), in which a fixed fraction f star ~ 1.7— 17% of the gas (normalised to reproduce the observed 

IGM metallicity at z = 3 of Z = 1CT 3 - W~ 2 Z Q ) in halos with T vir > 10 4 K (which are able to 

undergo atomic cooling) fragment to form a starburst lasting for ~ 10 7 years, and the halo collapse 

rate is given by the Press-Schechter formalism (Sasaki 1994): 

d fiform IdDdnps,^ , S 2 

where D{z) is the growth factor, and 5 C = 1.7 is the threshold above which mass fluctuations 
collapse. In this case the star formation rate is given by: 

1 1 Hi f'°° dNf° rm 

= — — ^ fstar / dt / dM ——(M,z)M (6) 

p c {z)t \l m Jt(z)-t J M(T vir =10*K,z) dM 

An approximate fit to the comoving star formation rate in the interval 3 < z < 10 can be given by 
A* = exp(a0 + alz + a2z 2 ) M yr" 1 Mpc~ 3 h 3 , where aO = -0.841, al = 0.395, a2 = -0.0295. The 
comoving emissivity is: 

e„(t) = p c [ dt'F u (t-t')n*(t') (7) 
J o 

where F l/ (At) is the stellar population spectrum, defined as the power radiated per unit frequency 
per unit initial mass by a generation of stars with age At. I obtain this spectrum from the Bruzual 
& Chariot (1993) code for a Z = 10~ 2 Zq population (i.e., extremely low metallicity), assuming 
a Salpeter IMF with lower and upper mass cutoffs at 0.1 and 100 Mq, and only computing the 
emissivity redward of the Lyman break. I ignore the effects of dust extinction, which should be 
negligible at these high redshifts when the metallicities are very low. 



3.1. Recombination radiation 

Is the number density of recombination line photons such as Lya sufficiently high to cause 
significant pair production opacity? It has been emphasized (e.g., Loeb k, Rybicki 1999) that apart 
from possible dust attenuation, Lya photons are not absorbed by the IGM but resonantly scattered 
until they redshift out of resonance. Indeed, the Lya radiation intensity has been predicted to be 
particularly strong prior to the epoch of reionization (Haiman, Rees &; Loeb 1997, Baltz, Gnedin 
& Silk 1998). This is because the optical depth in the Lyman series is very high, and all the 
Lyman series lines except Lya are absorbed immediately and redistributed. Thus, the Lya and 
Balmer lines are considerably brighter because they receive the energy of the Lyman series. This 
has spawned suggestions of detecting the epoch of reionization by a sharp drop in intensity at the 
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rest-frame Lya wavelength (Baltz, Gnedin & Silk 1998, Shaver et al 1999). In particular, Baltz, 
Gnedin Sz Silk (1998) find that Lya is about 3 times brighter and Ha is about 30 times brighter 
immediately prior to the reionization redshift. 

If this is indeed the case, Lya photons might contribute significantly to the pair production 
opacity. Let us define fj U mp = n recom i, t Lya/ ' n tot,Lyai the jump in the number density of photons 
longward of the Lya wavelength. This is given by the number of photons n recom b,Lya injected 
at the Lya wavelength , over the total number of photons n to t t Lya, include those redshifting into 
resonance. Most studies calculate the number of Lya recombination photons by summing over 
the IGM recombination rate in numerical simulations, by direct estimation of gas clumping in 
the simulations. In fact, this underestimates the number of Lya photons produced: if the escape 
fraction of ionizing photons is small, most recombinations occur in the HII regions of dense halos 
where star formation takes place, and the primary source of Lya photons are the ionising sources 
themselves. If each ionizing photon is converted into a Lya photon (plus lower energy photons) at 
the source, then h^ya ~ «ion(l — fesc), where f esc is the average escape fraction of ionizing photons 
from a souce. Thus, inserting e(z^, z) = EL ya hi on {l — f e sc)S(v — VLya) into equation (4), I obtain 
for v < VLya, the solution 

Jv m {z) = ^7rM — J"' 7— ! 3 ^ nW2 ^)(l - fesc) (8) 

4vr H u Lya (tt m (l + z s ) s + Q A ) 1/2 

where rii on (z) is the production rate of ionizing photons in proper coordinates, prior to attenuation 
by the ionizing photon escape fraction. This is a lower bound on J„ ya (z) since it does not include the 
conversion of photons trapped in higher order Lyman resonances to Lya photons. I find that for the 
adopted stellar spectra, in which the Lyman edge is typically f break ~ 5, then typically fj ump ~ 1, 
which implies that Lya photons are comparable to UV continuum photons as a source of opacity. 
The jump factor fj urnp may easily be understood as the ratio of the total number of Lyman alpha 
photons produced (or, the total number of ionizing photons produced) against the total number 
of photons in the 10.2 — 13.6eV range. If the break at the Lyman edge is reduced, then fjump is 
increased and Lya photons are more important in contributing to the pair-production opacity. In 
particular, for a power law spectrum with no discontinuity at the Lyman edge (e.g. as for quasars), 
Lya photons are the dominant source of opacity. Lya photons are also the dominant source of 
opacity in scenarios where the universe is reionized by zero metallicity stars. The higher effective 
temperature of these Pop III stars imply that their spectrum is much harder, and significantly fewer 
photons are emitted longward of the Lyman break (Tumlinson h Shull 2000). 

Lya photons could thus an important source of pair production opacity provided the escape 
fraction is small, i.e. (1 — f esc ) « 1, and dust attenuation while the Lya photons resonantly scatter 
in the IGM is unimportant. If Lya photons are the dominant source of pair production opacity, 
this would be extremely interesting, as it would provide a indirect census of the ionizing photon 
emissivity, prior to attenuation within the host sources. This could prove to be a good measure of 
the comoving star formation rate. One way to check if this is the case would be to directly image 
sources in NGST in rest frame UV longward of Lya (sufficiently far away from the Lya damping 
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wing), as well as in rest frame Balmer line emission (which should also be directly proportional 
to the production rate of ionizing photons, hua c< ^Lya °c hi on (l — f esc )). This will give a sense 
as to whether Lya or UV continuum photons are a greater source of opacity, provided the trend 
observed in bright sources extrapolates down to lower luminosities. 

3.2. Results 

Fig 2 shows the predicted attenuation factor as a function of observed photon energy, for 
grazars at redshifts z s = 3,6, 10, for both high (f s t ar = 17%) and low (f s t ar = 1-7%) star formation 
efficiencies. There are several important features to note. Firstly, the shape of the attenuation 
curve for Lya is similar to that for UV continuum photons. It is not possible to immediately 
distinguish between scenarios where Lya photons and UV continuum photons are the dominant 
source of opacity. Secondly, for lower star formation efficiencies, attenuation both sets in at only at 
higher energies and the attenuation curve is significantly shallower (i.e., it reaches full attenuation 
after a much longer energy interval). This is easy to understand. Gamma-ray photons of higher 
energy have both a longer path-length to travel before they redshift to E < E t h, and a higher 
number density n(eth < e < ei) of photons to pair produce against, since e t h ~ (m e c 2 ) 2 / 'E is 
lower. Thus, if the overall number density of UV photons is lower, one must go to higher energies 
to achieve the same attenuation. In Figure 3, I show the relative contribution of different rest- 
frame target photon energies and redshift intervals for gamma-ray photons with t(E ,z s ) = 1, for 
a variety of source redshifts. The shape of these curves is largely dependent on the overall number 
density of target UV photons. For a lower SFR, the energy E at which t(E ,z s ) = 1 increases, 
and the target photon interval and redshift interval contributing to the resultant opacity broadens. 
Note that most of the opacity comes from redshifts comparable to that of the source. Also, the 
opacity arises from a fairly narrow target photon energy interval ~ 3 — 13.6eV; this varies weakly 
with source redshift. Thus the results do not depend strongly on the assumed spectral slope of the 
UV sources. I also obtained by direct computation the contribution to the opacity due to ionizing 
photons shortward of the Lyman edge. Depending on one's assumptions for fesc, f break, /opacity, it 
is smaller by 2-4 orders of magnitude, and is completely negligible. 

An important technique for isolating the contribution of high-redshift radiation fields would 
be to use measurements of differential absorption between blazars at different redshift. Consider 
two blazars with at redshifts z\ and Z2, with Z2 > z\. At a given observed energy E , photons 
from each blazar will encounter identical optical depths for z < Z\. Thus, any additional opacity 
seem in the spectrum of the blazar at Z2 must arise from the redshift interval z\ < z < Z2 alone, 
At(E d ) = dzf^ j\d(cos9){l - cosO) J^den^e, z)a(E,e,6). Note that this identification is 
independent of any uncertainties in the spectral slope or redshift evolution of the UV background. 
As argued above, radiation field intensity fluctuations are unimportant since the pair-production 
opacity arises on scales of order a Hubble length. Thus, in Figure (2), any difference between the 
z s = 3, 6, 10 curves is solely due to radiation fields between 3 < z < 6 and 6 < z < 10; if star 
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formation ceases for z > 3, the curves would lie on top of one another. If there is little difference 
between two absorption curves from sources at z\ and Z2, one can get an upper bound on the UV 
emissivity in the range z\ < z < Z2', likewise, if two curves are so widely separated that meaningful 
measurements of At(E q ) cannot be obtained (in particular, if absorption saturates in one of the 
curves), one can get a lower bound on the UV emissivity in the range z\ < z < z 2 - In Figure (4), 
I show how the attenuation at a fixed observed photon energy E D is expected to rise with redshift, 
due to the opacity provided by high redshift UV fields. For lower star formation efficiency, the 
attenuation rises more slowly; if there is no star formation at high redshift the curve would be flat. 
In Figure (5) I show how the energy E for which t(E , z s ) = 1 is expected to fall with increasing 
redshift; due to the increased opacity provided by high redshift photons, attenuation sets in at 
lower energies. Again, if there were no star formation at high redshift, this curve would show no 
evolution. 

In fact, differential absorption measurements provide a powerful cross-check on the technique: if 
t(E q , Z2) < t(E , zi), then the assumed unabsorbed blazar spectrum must be evolving with redshift, 
due to a change in internal absorption, or underlying spectral index. Note that if star formation 
rates are high, the attenuation occurs rapidly over a small energy interval, and uncertainties due to 
the extrapolation from the unabsorption portion of the spectrum is small; for lower star formation 
rates the attenuation occurs over a larger energy interval, and uncertainties due to extrapolation 
are greater. 

4. Conclusions 

In this paper, I have suggested that if blazars can be detected at high redshift, detection 
of gamma-ray absorption due to pair production against high-redshift UV photons will provide a 
valuable probe of high-redshift UV radiation fields. This is because the sharp Lyman edge in the 
intergalactic radiation fields implies that gamma-ray photons have only a limited redshift interval in 
which to pair-produce. As they redshift to lower energies, they require photons with e > 13.6eV to 
pair-produce, so the universe becomes optically thin. The shape of the attenuation curve is primarily 
sensitive to the overall number density of photons longward of the Lyman edge at high redshifts: 
the higher this number density, the lower the gamma-ray photon energy at which pair-production 
opacity sets in. This makes it a useful test of the overall level of star formation and ambient UV 
radiation fields present at high redshift. Lya photons provide an important contribution to this pair 
production opacity, and indeed may be the dominant source of opacity if sources with a relatively 
lower fluxes longward of the Lyman edge (such as quasars or low metallicity stars) are abundant. 
Finally, measurements of differential absorption between blazars at the same observed energies will 
allow us to cleanly isolate the increase in opacity due to radiation fields at high redshift. 

There are two large uncertainties. The first is whether GLAST will be able to see high-redshift 
blazars at all. However, in the unified model of AGN, the scarcity of ultra-luminous blazars is a 
geometrical effect (due to relativistic beaming) rather than a requirement of extremely large black 
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hole masses. In fact, the luminosity boost provided by beaming reduces the black hole mass by 
several orders of magnitude below that demanded by the Eddington limit. So it is at least plausible 
that high redshift blazars will be detectable. The second uncertainty is whether absorption seen in a 
blazar will be internal, rather than due to pair production against photons in the IGM. However, at 
GeV energies we have some physical understanding of the observed EGRET spectra (e.g., Ghisellini 
et al 1998); opacity to gamma-ray photons due to internal radiation fields can be constrained by time 
variability arguments and other constraints. Furthermore, GLAST should assemble an extremely 
large catalog (> few thousand) of low redshift blazars, whose spectra can be studied in detail (and 
the contribution to opacity due to low redshift star formation can be quantified by other means); 
provided blazar properties do not evolve too strongly with redshift, we should have a firm handle 
on the intrinsic unabsorbed blazar spectrum. 
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Fig. 1. — Top panel: Best fit model(solid line) for blazar detection over the entire sky with EGRET 
detection threshold S(E > lOOMeV) = 2 x 10~ 7 photons s _1 cm~ 2 , against actual number of sources 
detected (dotted line). The model assumes 7 = 6 and L™$ nsic = 3.2 x KT 2 L^ 9 13 . Bottom 
panel: predicted number of blazar detections with GLAST detection threshold S(E > lOOMeV) = 
2 x 10~ 9 photons s _1 cm -2 as a function of redshift. 



UV continuum 




Fig. 2. — The predicted attenuation factor as a function of observed photon energy E Q for UV 
continuum opacity (top panel) and Lya photon opacity only, for source redshifts z s = 3, 6, 10. The 
solid lines are for a high star formation efficiency (f s tar = 17%), and dashed lines are for a low star 
formation efficiency (f s tar = 1.7%). 
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Fig. 3. — The contribution of rest-frame target photon energy interval (top panel) and redshift 
interval (bottom panel) to the opacity, all for gamma-ray photons for which t(E , z s ) = 1. Figures 
are for 3 different source redshifts (z s = 3, 6, 10) and 3 different UV emissivity scenarios:high SFR 
(fstar = 17%, dark solid line), low SFR (f s tar = 1-7%, dashed line), Lya photon opacity only 
(assuming high SFR, thin solid line). Note that contribution of different energy intervals evolves 
only weakly with source redshift (z s = 3,6,10 from right to left in top figure). As the number 
density of target photons decreases (low SFR), the distribution of energy and redshift interval 
contributions broadens. 
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Fig. 4. — The attenuation factor as a function of blazar redshift for a given observed photon 
energy. Solid curves are for the high star formation efficiency case, dashed curves are for the low 
star formation efficiency case. Any increase in the attenuation between z\ and Z2 is solely due to 
UV radiation in the range z\ < z < Z2- Note that if the star formation efficiency is low, the degree 
of attenuation changes much more slowly with redshift. 
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Fig. 5. — The variation of the observed energy E Q at which r = 1 with source redshift, for a 
high star formation efficiency {fstar = 17%), a low star formation efficiency (f s tar = 1-7%), and a 
high star formation efficiency model where Lya photons provide the main source of opacity. Note 
how the curve flattens at high redshift due to the reduced opacities at high redshift. If the pair 
production opacity at high redshift is negligible, the curve would be completely flat. 



